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Broadband suppression of quantum noise below the Standard Quantum Limit (SQL) becomes 
a top-priority problem for the future generation of large-scale terrestrial detectors of gravitational 
waves, as the interferometers of the Advanced LIGO project, predesigned to be quantum-noise- 
limited in the almost entire detection band, are phased in. To this end, among various proposed 
methods of quantum noise suppression or signal amplification, the most elaborate approach implies 
a so-called xylophone configuration of two Michelson interferometers, each optimised for its own 
frequency band, with a combined broadband sensitivity well below the SQL. Albeit ingenious, it is a 
rather costly solution. We demonstrate that changing the optical scheme to a Sagnac interferometer 
with weak detuned signal recycling and frequency dependent input squeezing can do almost as good 
a job, as the xylophone for significantly lower spend. We also show that the Sagnac interferometer is 
more robust to optical loss in filter cavity, used for frequency dependent squeezed vacuum injection, 
than an analogous Michelson interferometer, thereby reducing building cost even more. 


I. INTRODUCTION 

The past decade was marked by great achievements 
in gravitational-wave (GW) instrumental science. The 
international network of gravitational wave (GW) de¬ 
tectors, comprising three LIGO (Laser Interferometer 
Gravitational-wave Observatory) detectors [1-3] and the 
EGO (European Gravitational Observatory) detector 
Virgo [4, 5] has reached the project sensitivity and has 
accomplished 6 runs of scientific data collection. Apart 
from setting limits on a population of various sources of 
gravitational radiation in our Galaxy and beyond [6-11], 
all detectors in the network have reached hitherto un¬ 
seen displacement sensitivity of ~ 10“^® m/-\/Hz, making 
them, perhaps, the most sensitive displacement sensors 
in the world. The figure 10“^® m/-\/Hz is remarkable by 
itself as it is only 1 order of magnitude above the magni¬ 
tude of quantum zero-point fluctuations for a mechanical 
object with a mass of 10 kilogramms that is the mass of 
the LIGO interferometer core optics mirrors. 

The second generation detectors, like Advanced LIGO 
[12, 13], Advanced Virgo [14], KAGRA [15] and GEO-HF 
[16] will be already operating at the level of this quantum 
limit, meaning that their sensitivities will be governed 
by quantum fluctuations of light in the frequency band 
around 100 Hz wherein most of the target GW sources 
are expected to emit. Detection band below 100 Hz will 
be dominated by quantum fluctuations of light ampli¬ 
tude, known as quantum radiation pressure noise (RPN), 
while at higher frequencies detector sensitivity will be 
limited by quantum phase fluctuations, usually referred 
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to as quantum shot noise (SN). The best sensitivity point, 
where these two noise sources become equal, is known as 
the Standard Quantum Limit (SQL) [17], which, in a 
broader context, characterises the regime in which the 
quantum measurement noise (SN) becomes equal to the 
back action noise (RPN) — the latter one being a direct 
consequence of the Heisenberg uncertainty principle. 

The design and configuration of the third generation 
detectors remain an open question. However, there is a 
consensus in the community that they should have an 
order of magnitude better quantum noise than their pre¬ 
decessors [19, 20] in as broad frequency band as possi¬ 
ble. At the same time, this sensitivity gain should be 
achieved in as economic way as possible, which implies 
that a single-interferometer solution should be favoured 
over any multi-interferometer ones. In this article, we ar¬ 
gue that the most efficient way towards broadband quan¬ 
tum noise suppression goes via the change of interferom¬ 
eter configuration. We show that Sagnac interferome¬ 
ter has superior potential for broadband sensitivity gain 
compared to Michelson interferometer for any given set of 
advanced interferometric techniques, i.e. signal recycling 
[21-23], squeezed vacuum injection [46, 47], frequency- 
dependent phase rotation of squeezed vacuum phase [24] 
etc. We also demonstrate that better broadband per¬ 
formance is exhibited by Sagnac interferometer with sig¬ 
nificantly relaxed requirements to auxiliary optics and 
thereby at much lower cost. 

In Fig. 1, we show quantum noise sensitivities for two 
variants, namely for a LIGO-sized interferometer with 4- 
kilometre arms and 40 kilogram mirrors, and for the in¬ 
terferometer with parameters planned for Einstein Tele¬ 
scope, i.e. for an interferometer with 10-kilometre arms 
and mirrors of 200 kilogram. Full list of optical parame¬ 
ters for these configurations that was obtained by means 
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FIG. 1: Optimal quantum noise curves of Sagnac (red line) and Michelson (blue dashed line) interferometers are shown for 
two variants, i.e. LIGO-like configuration {Left panel) features arm length L = 4 km, mirror mass, M = 40 kg, and maximal 
circulating power in the arm, = 1 MW, while ET-like configuration {Right panel) features arm length L = 10 km, mirror 

mass, M = 200 kg, and maximal circulating power in the arm, = 3 MW. In both cases, we consider a set of QND 

techniques that are currently considered as the most plausible ones for implementation in the future detectors, i.e. frequency- 
depended squeezed vacuum injection (filter cavity loss Af/L[ = 1 ppm/m), balanced homodyne phase readout and detuned 
signal recycling. Optimisation of interferometer parameters is done so, as to match closely a broad-band lower bound quantum 
noise curve (solid black line) that is the best possible quantum noise sensitivity one can get in principle for given level of optical 
loss, given circulating power in the arms and fixed interferometer bandwidth (see details and discussion in Sec. H. For scale, 
we show baseline design sensitivity curves of wide-band aLIGO {Left panel, dark green dashed line) [18] and of a xylophone 
configuration of the Einstein Telescope {Right panel, dark green dashed line) [19]. Parameters of respective optical setups are 
given in Table I. 


of optimisation procedure outlined in Sec. IV is given in 
Table I. These results answer the question: what is the 
best single interferometer configuration that has quan¬ 
tum noise as broadband and as low as the total noise of 
the aLIGO broadband baseline configuration [18] and of 
the ET-D xylophone configuration [19], respectively, for 
a given set of advanced interferometric techniques a.k.a. 
QND techniques, i. e. frequency dependent 20 dB squeez¬ 
ing, detuned signal recycling and balanced homodyne de¬ 
tection? Table I clearly demonstrates Sagnac scheme 
advantage over Michelson in an entire frequency band. 
Noteworthy is the fact that this excellent result can be 
achieved with a single filter cavity with rather realistic 
optical loss, 1 ppm/metre, that can be already achieved 
in the laboratory [25]. 


II. ON COMPARISON OF DIFFERENT 
INTERFEROMETER CONFIGURATIONS AND 
CHOOSING THE BEST ONE. 

Before we move on to details of quantum noise calcu¬ 
lations of considered interferometer schemes, it is impor¬ 
tant to have an agreement on how to rate those config¬ 
urations against their ability to detect GWs. Quantum 
noise of the 3rd generation interferometers is a compli¬ 
cated interplay of several advanced quantum techniques, 
mentioned above. It results in a quite complex depen¬ 
dence of the interferometer sensitivity on a multitude of 
parameters. This raises a question of optimisation of 
quantum noise curve and finding the best combination 
of those parameters. But what “the best” means in this 
context. How do we define the criterion that resolves 
what configuration is better than the others? 

Many different answers were given to this question in 
the literature. Some optimised the interferometer signal- 
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Parameter 

Notation 

Michelson 

Sagnac 

Mirror mass, kg 

m 

40 

200 

40 

200 

Arm length, km 

L 

4 

10 

4 

10 

Circ. power, MW 

Pc 

1.0 

3.0 

1.0 

2.7 

ITM transmittance, % 

Titm 

15 

5 

7 

11 

Squeezing factor, dB 

CdB 

13.2 

13.3 

20.0 

19.4 

(w/5% inj. loss) 


(10.2) 

(10.3) 

(12.3) 

(12.1) 

Squeezing angle, deg. 

0sqz 

-6° 

-4° 

-27° 

-24° 

Homodyne angle , deg. 

C 

84° 

86° 

63° 

67° 


Signal recycling cavity parameters 


SRM transmittance, % 
SRC detuning, deg. 

PSRM 

4‘SRC 

82 

90° 

12 

90° 

77 

102° 

67 

99° 

Filter cavity parameters 

FC detuning, Hz 

Sf 

34 

13 

525 

456 

rFC bandwidth, Hz 

7/1 

34 

17 

922 

767 

LFC mirr. trans., ppm/m 

Tf/Lf 

2.8 

1.5 

77.3 

64.3 

rFC loss bandwidth, Hz 

7/2 

12 

12 

12 

12 

LFC r. trip loss, ppm/m 


1.00 

1.00 

1.00 

1.00 


TABLE I: Optimal parameters for configurations of Sagnac 
and Michelson interferometers which quantum noise sensitiv¬ 
ity curves are presented in Fig. 1. 


to-noise ratio for specific GW sources [26, 27]. For in¬ 
stance, the most popular figure of merit in GW com¬ 
munity, the detection range, is nothing more than a 
renormalised SNR for the detection of GWs emitted by 
an etalon compact binary system (comprised of neutron 
stars or black holes) in the course of inspiralling phase of 
their evolution. More general criterion was proposed in 
[28] , where the optimisation seeks to provide the broadest 
possible total noise curve (a sum of a quantum noise and 
other classical noise sources), thereby seeking to include 
as many various astrophysical sources in the detection 
band of the antenna, as possible. In both cases, how¬ 
ever, the optimisation depended heavily on the model of 
non-quantum noise for a specific project, i.e. aLIGO, or 
Einstein Telescope. 

We argue that this approach is very restrictive and 
does not allow to reveal the potential of different inter¬ 
ferometer configurations in full, for it is squeezed into 
a narrow frame of existing noise models. Nevertheless, 
we know examples when a new technology comes into 
play and predictions of classical noise models have to be 
revised dramatically [29]. 

In this paper, we suggest that quantum noise of dif¬ 
ferent configurations shall be compared against the lim¬ 
itations of quantum origin. Namely, we suggest that a 
lower bound of quantum noise, which all considered con¬ 
figurations shall be compared against, is to be derived as 
a sum of limitations imposed by two main parameters for 
any quantum noise-limited interferometer, i.e. the level 
of optical loss in the main optics and the finite optical 
power circulating therein. 


Incoherent loss-associated vacuum fields that enter 
the interferometer in accordance with the fluctuation- 
dissipation theorem [30] create random radiation pres¬ 
sure force that cannot be compensated by any quantum 
technique. If, due to various reasons, the total fraction 
eioss of all photons that enter the interferometer is lost 
therein, and if the effective squeezing of vacuum injected 
in the dark port is given by a factor [see eq. (47)], 

the ultimate limit for residual radiation pressure noise (in 
GW strain spectral density units) reads (cf. Eq. (413) in 
[31]): 


‘S'rp loss(^) — ^SQL(^)£lois 6 ) ( 1 ) 

where 

^SQL(f^) = y^^||^, (2) 

is the Standard Quantum Limit of the interferometer in 
terms of GW strain^ We will use this limit as a first 
component of our QN lower bound curve. 

High frequency region of all interferometers is domi¬ 
nated by shot noise, or quantum phase fluctuations. Its 
rise on upper frequencies is determined by a finite band¬ 
width, 7 , of the detector. For a simple resonance-tuned 
Fabry-Perot interferometer with squeezed vacuum and 
lossy optics, a shot noise contribution can be written as: 


s!^ 


SN 


(n) = % 




■ ^loss 


/Ceff(f^) 


( 3 ) 


where /Ceff is the frequency-dependent optomechanical 
coupling strength introduced by Kimble et al. [24] and 
equal to: 


^eff(^^) = 


0T 


2T'eff 


1 -I- 2\/Rcs cos 2Qt + i?eff 


( 4 ) 


with 0 = AlooPcI{McL) and Pc is the total optical power 
circulating in the interferometer, t = L/c is light travel 
time between the mirrors of the cavity and Poff = 1 ~ Feff 
is the effective reflectivity of the cavity. According to 
scaling law derived by Buonanno and Chen in [32], the 
Fabry-Perot interferometer is fully equivalent to a signal- 
recycled Fabry-Perot-Michelson interferometer in terms 
of quantum noise with (in the the only substitution of 
Pc = Peff where Peg stands for light power circulating in 
the first one. 

For fixed bandwidth 7 = cTes/4L = 2 tt ■ 1000 s“^, cir¬ 
culating power Pc and achievable effective squeezing at 
the dark port of the interferometer, Eq. (3) sets 


^ Here we used a particular formula for the SQL of an interfer¬ 
ometer with 4 test masses of equal value M and the mechanical 
dARM-mode defined as x^arm = (a:BTM “ ^itm) “ (^ETM “ 
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the ultimate high frequency limit on quantum noise for 
interferometers, not using active optomechanical ampli¬ 
fication techniques [33]. Therefore we will use it as the 
second component of our QN lower bound curve, which 
can now be expressed in terms of lower bound quantum 
noise spectral density as: 

‘S'ref(^) = '5'rp loss(^) + 5 'sn(^) • (5) 

The resulting lower bound quantum noise curves (for 
LIGO-like and for ET-like interferometers) are plotted in 
Fig. 1 and will be used throughout the rest of the paper as 
yardsticks for different configurations of interferometers. 

III. QUANTUM NOISE OF THE 
SIGNAL-RECYCLED SAGNAC 
INTERFEROMETER 

In this section, we concentrate solely on quantum fluc¬ 
tuations of light and the influence they have on the sensi¬ 
tivity of the GW detector based on Sagnac interferometer 
topology. 


double-measurement variant of the quantum speed me¬ 
ter, shown in Fig. 3. Indeed, visiting consequently both 
arms (see Fig. 2), counter propagating light beams ac¬ 
quire phase shifts proportional to a sum of arms length 
variations XN,E{t) = (a^E™(^) ~ 2 :^^(t)) (hereinafter 
I(E)TM stands for Input (End) Test Mass) for of both 
cavities taken with time delay equal to average single 
cavity storage time Tarm: 

S(j)R OC XNit)+XE{t + 

'^arm) ; (7) 

5(pL OC XEit) + XN{t + '^arm )• (8) 

After recombining at the beam splitter and photo detec¬ 
tion the output signal will be proportional to the phase 
difference of clockwise (R) and counter clockwise (L) 
propagating light beams: 

SipR - 54 >l OC [a;7v(t) - XNit + Tarm)]- 

- [XEit) - XE{t -I- Tarm)] OC 

OC XAr(t) - i;E(t) + O(rarm) (9) 

that, for frequencies fl <C Ta^m; is proportional to relative 
rate of the interferometer arms length variation. 


A. Suppression of radiation pressure noise in 
Sagnac interferometer. 


B. Transfer matrix for calculation of quantum 
noise. 


The power of Sagnac interferometer (SI) to reduce 
back-action noise is nested in its ability to sense the 
relative rate, or in other words speed of an arm cavi¬ 
ties length variation, whereas Michelson interferometer 
senses arms length variation itself. The simple way to 
understand how a speed measurement can reduce back- 
action is to consider a simple thought experiment de¬ 
picted in Fig. 3. Here the free mirror is sensed twice 
by the same laser light that is reflected from both the 
front and the rear surfaces thereof with a time delay r 
between reflections. The phase of outgoing light is mea¬ 
sured by, say homodyne detector, and is proportional 
to the the difference of the succesive mirror coordinates: 

OC (x(t -I- r) — x{t)) ~ ur, where v stands for the 
mean velocity of the mirror over the interval r. If the 
signal force one seeks to measure, watching the change of 
the mirror velocity, has characteristic frequency Q much 
smaller than the two kicks light gives to the mir¬ 
ror on the consecutive reflection partly compensate each 
other and the resulting back-action force turns out to be 
depressed by a factor oc Hr <C I: 




—iflr 


2Pr 


pulse 


( 6 ) 


As shown by Gaves and Schumaker in [40, 41], quan¬ 
tum noise of light in any parametric optical device can be 
conveniently described within the frames of two-photon 
formalism. Namely, noise can be considered as tiny 
stochastic variations of quadratures of carrier field travel¬ 
ling through the device. Any variations of interferoeme- 
ter parameters due to signal force, e.g. differential arms 
length change in GW interferometer, also lead to varia¬ 
tion of outgoing light quadratures and thus can be de¬ 
scribed using the same formalism. 

In two-photon formalism, one starts with writing the 
input and outgoing light fields of the interferometer at 
some fixed location via its sine and cosine quadratures: 


F“(t) 

F°“‘(t) 


= £0 [(^™ + a™) cos Wpt-I-a*" sinwpt] , (10) 
= £0 (i?°“‘-I-S°“‘) cos Wpt-I-sinwpt ().l) 


where A*” (i3°“*) is classical mean amplitude of the in¬ 
put (output) light at pump laser frequency ujp and a™^ 
(6™ ) are small quantum variations and variations due 
to signal force that can be defined through creation and 
annihilation operators as 


as compared to the back-action of single light pulse with 
an average power Ppuise which one expects in a single 
reflection experiment sensitive to the test mass displace¬ 
ment. 

As was realised by Ghen and Khalili [34, 37], that 
the zero area SI [38, 39] actually implements the initial 


he 


a -I- 


and Os 


a — 

iV2 


( 12 ) 


and similarly for outgoing fields. Note that we do not 
specify time as an argument here, as the same definition 
holds in spectral domain that we make use of in the rest 
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FIG. 2: Two possible optical realisations of zero-area Sagnac speed meter. Upper left panel: Sagnac interferometer with 
triangular ring cavities in the arms as proposed in the original work by Chen [34]. Upper right panel: Sagnac interferometer 
with linear Fabry-Perot cavities in the arms and polarised optics, consisting of additional polarisation beam splitter and A/4- 
plates, used to direct light from one arm to another as per proposal by Danilishin [35] and later analysed by Wang et al. in 
[36]. Lower panels show simplified schematics of frequency-independent {lower left panel) and frequency-dependent squeezed 
vacuum injection (lower right panel) and serve to introduce notations for input and output fields in both situations. 


f (!(?!) oc x{t -b t) — x{t) ~ v^T J 



FIG. 3: Thought experiment on test object speed measure¬ 
ment: The free mirror is sensed twice by the same laser light 
that is reflected from both the front and the rear surfaces 
thereof with a time delay r between reflections. The phase 
of outgoing light is measured by, say homodyne detector, and 
is proportional to the the difference of the successive mirror 
coordinates: (j>out oc {x{t + r) — x{t)) ~ vt, where v stands for 
the mean velocity of the mirror over the interval t. 


Fourier transform as follows: 

-oo 

So, in order to fully describe signal and noise in a (loss¬ 
less) GW interferometer it is sufficient to track the trans¬ 
formation it does to the quadrature operators of input 
light that can be represented in a matrix form as trans¬ 
formations of 2-dimensional vectors a = {dc, dg}"'" and 
b = {Sc, Ss}"'" and GW signal h{n). For GW detector it 
can be written in general form as: 

S = T ■ d -I- t , (13) 

^SQL 

where 

[Tcc(O) Tc«(F!) 

[T,c(f^) T,s{^)_ 

is the optical transfer matrix of the interferometer, 

tc{fl) 

ts{n) 


of this article. These two domains are connected through 


( 15 ) 
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is an optical response of the interferometer on an ex¬ 
ternal influence that is denoted as and d:sQL is the 
corresponding SQL for the mechanical degree of freedom 
of the interferometer expressed in units of X. In GW in¬ 
terferometry X is either the displacement of the mirrors, 
X, or the external force, F, that causes this displacement, 
or, more habitually, the GW strain, h. In each case, the 
corresponding SQL applies. The relation between these 
three quantities is discussed in Sec. 4.3 of [31]. 

Note, that the I/O-relation form (13) means for optical 
detuning-less system the following expression: 


6 = 


1 0 

-jc{n) 1 


a -f e 


i/3(n) 


0 


X 

-^SQL ’ 


where /C is a Kimble optomechanical coupling factor and 

is some frequency-dependent phase shift optical fields 
acquire when passing through the interferometer and 
that does not affect the expression for spectral density 
of quantum noise. 

If one assumes the interferometer readout quantity oq 
(proportional to homodyne detector photo-current 
as a result of a homodyne detection with a homodyne 
angle C, i.e: 


= be cos C + ds sin (^ = HJ • b , = 


cosC 

sin() 


, (16) 


the spectral density of quantum noise at the output port 
of the interferometer can be obtained using the following 
simple rule: 


S^i^) — dsQL 


HT . T • S“ • Tt • 

IHj-t.P 


(17) 


where S“ stands for spectral density matrix of injected 
light and components thereof can be defined as: 


i(m|di(II)(dj(II'))l -I- (aj(il'))lai(r2)|In) , (18) 

where \in) is the quantum state of vacuum injected in 
the dark port of the interferometer and (i,j) = {c, s} 
(see Sec. 3.3 in [31] for more details). In present article 
we deal with single-sided spectral densities S and hence 
in case of input vacuum state: 

\in) = |uac) ^ S“ = ^ ® . 

Now with all necessary tools in hand we can proceed 
with the analysis of Sagnac interferometer and derive 
spectral density of quantum noise for it. 


C. Sagnac interferometer I/O relations. 


definiteness, in this section, we stick to a configuration 
of Sagnac interferometer that utilises ring arm cavities 
(as per the left panel of Fig. 2), although the results we 
obtain are applicable to both realisations unless loss is 
taken into account. 

Unlike Michelson interferometer, in Sagnac interferom¬ 
eter light beam visits two arm cavities before recombina¬ 
tion with a counter-rotating beam at the beam splitter 
(see Fig. 2). At the same time, two light beams hit the 
cavity, one coming directly from the beam splitter and 
the one, that has just left another arm. In notations of 
Chen’s paper [34] quadrature operators of light enter ing 
and leaving the arm can be identified with two indices JJ, 
e.g. , where I stands for the either of two beams, L or 
R, and J stands for the either of two arms (J = E,N). 
Here R marks the light beam that first enters North arm 
and then travels the interferometer in the right direc¬ 
tion (clockwise), and L marks the beam travelling the 
interferometer in the opposite (counterclockwise) direc¬ 
tion after entering the interferometer through the East 
arm. Thus, single lossless arm I/O relations read, as¬ 
suming high-finesse arm cavities (Tixm ^ for general 
case see Appendix C): 

bU _ g2z/3„„(n)^7J ^ 

(^)[af-/Carm(a" + al'^)] 

( 20 ) 

«SQlL 


with / indicating the other beam than I, i.e. R = L and 
L = R, hj = X— Xj”"^ is the arm elongation induced 
by signal force {e.g. gravitational wave tidal force), and 


Aa.rm — 


0a 


2Tj 


ITM 


1 — 2y/ i?iTM COS 2Ur -|- i?iTM 

20a 


•'arm'Tarm. 


+ H2) ’ 


/3arm = arctan 


1 -l- Vi?iT: 


'M 


1 — \/i?lTM 

~ arctan (U/7arm), 


tan Ht I ~ 


( 21 ) 


( 22 ) 


with 0arm = dwoFarm/(AfcL) and Farm = Pc/4:, where 
Pc is the total optical power circulating in both arms 
and 7arm = T'itm/(4t) is the half-bandwidth of an arm 
cavity. The final, approximate expressions above are 
obtained assuming that cavity linewidth and signal fre¬ 
quency are much smaller than cavity free spectral range 
^FSR = (27")This approximation nearly breaks down 
for detectors with arm length > 10 km, like Einstein Tele¬ 
scope, at frequencies of the order of 10 kHz, therefore we 
present exact formulae as well. 

Then, it is straightforward to derive full Sagnac I/O- 
relations, using junction equations for the fields at the 
output beam splitter (ring-cavity topology): 


First we consider a bare lossless zero-area Sagnac inter¬ 
ferometer and derive its input-output (I/O) relations.For 


-RN 


p + i 


-LE 

a 


p — i 


0 = 


-LN -RE 

b -b 

V2 


, (23) 
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as well as continuity relations between the beams that 
leave one arm and enter the other: 


down the I/O-relations for the lossless signal recycled in¬ 
terferometer in matrix form: 


'RE 


= b 


RN 


'LN 


= b 


'LE 


(24) 


The resulting I/O-relations for lossless zero-area Sagnac 
interferometer then read: 


be 

— g2i/3si 

1 

0 


A 

bs 


1 

1 

1 


is 


0 

V2/Csi 


i/3si ^ 

^SQL ’ 

(25) 


with coupling constant /Csi defined as: 

JCsi = 4 /Carm sin" /?arm=:^ 

+ 7arm)^ 

where ©si = 40arm and additional phase shift: 

/3SI = + I (27) 

One can now calculate spectral density of quantum 
noise of the zero-area Sagnac, using Eq. (17), where 
transfer matrix T and response vector t read: 


T = 

1 0 

, t = 

0 


-ICsi 1 


•s/2/Csi 


(28) 


Therefore one gets this simple expression for spectral 
density (it is the same for all tuned interferometers with 
balanced homodyne readout of quadrature 6^ and vac¬ 
uum state at the dark port, save to the expression for 

icy. 


,h _ ^SQL j [/Csi - COtC]^ + 1 
“ 2 I /Csi 


(29) 


D. Signal recycling in Sagnac interferometer 


<PSR 



f 


1 
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■i 

6 

V_ 

_ J 


FIG. 4: Schematic of I/O-relations for a signal-recycled inter¬ 
ferometer. 

Consider now how an addition of signal recycling mir¬ 
ror in the dark port of the Sagnac interferometer changes 
its quantum noise. A scheme of signal recycling cavity is 
shown in Fig. 4. Essentially, it acts as a feedback loop for 
the light leaving the main interferometer, and its sign is 
defined by the phase shift signal light acquires in a round 
trip across the signal recycling cavity (SRC). The feed¬ 
back features of the SRC become evident if one writes 


OSR = TTsRisR + tSR(^/^SQL) , (30) 

where 

Tsr = -PsrI + esrK,^sr(II “ PSR'irK2<^sR)~^™dsR 

(31a) 

tsR = 'rsRl^<^sR(I ~ PSRTrK2<^sR)~^t (31b) 

with psR and tsr are SR mirror amplitude reflectivity 
and transmissivity (sign notations are given in Fig. 4), I 
is identity matrix and 


— 


cos ip — sin ip 
sin ip cos ip 


(32) 


is rotation matrix responsible for propagation of light by 
a distance l^p resulting in a phase shift ip = ujqI^/c. 

These expressions are universal and allow to obtain 
transfer matrices for any (lossless) signal-recycled inter¬ 
ferometer characterised by a bare transfer matrix T and 
response vector t. A general case of arbitrary SRC de¬ 
tuning and lossy interferometer is considered in Appen¬ 
dices C 3 and C 4. Here we only give expressions for two 
special cases known as Resonant Sideband Extraction 
(RSE) [21, 42, 43] and Tuned Signal Recycling (TSR) 
[22, 23] that stand for the cases of SRC tuned in anti¬ 
resonance and in resonance, respectively. 

In terms of control theory, the first one corresponds 
to a negative feedback, as the light reflected off the SR 
mirror enters the interferometer with the opposite sign 
(2 (/sr = tt phase shift in a full round trip in the SRC). 
The second case gives a positive feedback (light returns to 
the interferometer in the same phase). Intermediate case 
is more complicated as it involves modification of the test 
masses dynamics by means of creating an optical spring 
[44, 45]. 

Substituting Eq. (28) into Eq. (31) with ^SR = 0 
(TSR) and 7r/2 (RSE), one arrives at I/O-relations of 
a form identical to that of Eq. (28), but with modified 
optomechanical coupling strength and overall phase shift 
that read: 


K. 


SI,SR 


'^SR^SI 


■^SR^SI 


/3: 


'SI,SR 


1 ± /9sRe^*^®ip 
1 ± PSR 


= arctan 


= —arccot 


1 T PSK 
1 T PSR 
1 ± PSK 


11 T Psrc^*^”” I 

tan ,5si ) = 


tan 2/3a 


(33) 


(34) 


where we denote TSR case by ”+” and RSE case by 
in accordance with the sign of feedback it creates. 

In approximation of fl <C i^fsr, the expressions for /Cgj 
transform as: 


j^± _ _ 40317 arm F± 

si.SR + 


( 35 ) 


































where F± = (1±psr)/(1TPSr) ■ This approximate form 
clearly displays the effect of signal recycling on quantum 
noise: it makes the overall bandwidth of the interferom¬ 
eter wider at the expense of lower peak sensitivity, when 
feedback is negative, i.e. for c/isr = 7r/2 (RSE), and nar¬ 
rower, but with improved peak sensitivity, for positive 
feedback when = 0 (TSR). Fig. 5 show the influence 
of signal recycling mirror on the quantum noise sensitiv¬ 
ity of both, Michelson (panel (a)) and Sagnac (panel (b)) 
interferometers. 


E. Ponderomotive squeezing and radiation 

pressure noise in recycled Sagnac interferometer. 

It is instructive to investigate how signal recycling ef¬ 
fects the frequency dependence of K. for Michelson and 
Sagnac interferometers and thereby radiation pressure 
noise in both interferometers. The appearance of K, in 
the off-diagonal component of the light field transfer ma¬ 
trix T is responsible for the transformation of the quan¬ 
tum state of light, passing through the interferometer, 
known as ponderomotive squeezing. Kimble et al. [24] 
showed that the transform (25) is equivalent to a se¬ 
quence of phase plane rotation followed by a squeeze and 
finally concluded by another rotation (see Appendix B 
for details). Vector of quadratures of the input field, 
{he, ds}, undergoes a clockwise rotation by a frequency- 
dependent angle t^pond = ?i>pondj then gets anti-squeezed 
(squeezed) along the Oc (os) quadrature axis by a factor 
gj-po„d (e“'’p°nd), and finally is rotated counter-clockwise 
by an angle Upond = '('pond + ^'pond, where 

fC 

^'pond = arctan — , (36) 

and 

e’^pond = -f ^ , <()pond = ^arccot^ . (37) 

Using a trigonometric identity, arctan a: -I- arccota: = 7r/2 
for a: > 0, it is straightforward to show that t'pond ~ 
Mpond = 7r/2. Frequency dependences of ponderomotive 
squeezing factor, rpond, and angle, Mpond, are given in two 
lower rows of Fig. 5 and illustrated by the transformation 
of the noise ellipses of the vacuum field at the input of 
the interferometer to the ponderomotively squeezed vac¬ 
uum at the output thereof (low frequency noise ellipse 
for Michelson case is drawn not to scale to fit it into a 
figure). 

One can readily see the feedback effect of the signal 
recycling cavity on the level of radiation pressure noise. 
Negative feedback (RSF) expectedly makes ponderomo¬ 
tive squeezing smaller in both, Michelson and Sagnac 
cases. Positive one (TSR) has an opposite effect. The 
superior radiation pressure noise performance of Sagnac 
interferometer is vivid from these plots too. 

Finally, we can see that squeezing angle, Mpond; varies 
less in Sagnac than in Michelson. This, as we will see 


later, allows to relax requirements for filter cavities for 
frequency dependent squeezing angle rotation in Sagnac 
case. 


F. Injection of squeezed vacuum and frequency 
dependent squeezing. 

Injection of squeezed vacuum into a signal port of the 
interferometer in order to reduce its quantum noise has 
been first proposed by Unruh [46] and then further de¬ 
veloped by Caves [47]. 

a. Constant squeezing. Squeezed vacuum used in 
this technique is a product of nonlinear interaction of 
light with a Kerr medium of a nonlinear crystal. The 
crystal is pumped at double carrier frequency 2a;p, and 
vacuum photons at carrier frequency are produced from 
the pumping photons in the process known as parametric 
down conversion (PDC). Since there is a strong correla¬ 
tion between the produced photons, quantum fluctua¬ 
tions of the vacuum field turned out to be ’’squeezed” in 
one quadrature and is amplified by the same factor (ide¬ 
ally) in an orthogonal one. From the mathematical point 
of view this transformation is similar to the ponderomo¬ 
tive squeezing and thus can be described by expressions 
given in Appendix B. For now, we are interested in the 
spectral density matrix, of a squeezed vacuum with 
arbitrary squeezing factor r and angle A. It is readily 
obtainable from Fq. B1 and taking into account that in 
a vacuum state \in) = \vac) the spectral density ma¬ 
trix (18) equals to identity matrix, one can calculate the 
spectral density matrix for squeezed light from (B3) as: 

= Ma • §[2r] • . (38) 

Substituting it into (17) and using Fqs. (28), one obtains 
a universal expression for quantum noise spectral density 
of an interferometer with injected squeezed vacuum: 

{e2’'(sin A - cos A [A - cot C])^+ 

-l-e“^’'(cos A-I-sinA [/C — cot^])^} • (39) 

This expression is valid for any interferometer with a 
transfer matrix and a response vector represented in the 
form of Fqs. (28), meaning for us, it is true for both, 
Michelson and Sagnac configurations, with or without 
resonance signal recycling. 

Injection of squeezed vacuum pursues a clear goal 
of reducing quantum noise in the readout quadrature. 
Were it not for radiation pressure noise and pondero¬ 
motive squeezing, the optimal squeezing angle would be 
A = C-|-7r/2, resulting in lower shot noise. For exam¬ 
ple, for phase quadrature readout, ^ = 7r/2, one requires 
a phase-squeezed vacuum, that is A = 0. 

Radiation pressure noise, however, makes constant 
squeezing not that beneficial. Imagine, we feed the 
Michelson interferometer with phase squeezed vacuum. 
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FIG. 5: Effect of signal recycling on Michelson (left) and Sagnac (right) interferometers quantum noise sensitivity (first row, 
panels (a) and (b)), ponderomotive squeezing strength (middle row, panels (c) and (d)) and quantum state rotation angle, 
rtpond (bottom row, panels (e) and (f)). Plots are drawn for ET typical parameters: arms length L = 10 km, mirrors mass, 
m = 200 kg, arm circulating power, Parm = 3 MW, phase quadrature readout, ^ = tt/2. ITM power transmissivity for Michelson 
is equal to Titm = 5.2%, and for Sagnac is Titm = 9.1%. Light quantum state transformation by an interferometer (w/o signal 
recycling) is illustrated by noise ellipses. Circles on panels (a) and (b) stand for injected vacuum noise ellipse. Ellipses on panels 
(e) and (f) demonstrate quantum state of the outgoing light transformed by ponderomotive squeezing at different frequencies 
(10, 100, and 1000 Hz). Green arrows point at the direction of readout quadrature, ^ = tt/2. 


Ponderomotive squeezing will transform it according 
with plots shown in panels (c) and (e) of Fig. 6. At 
higher frequencies (where /C is close to 0) there are no 
changes to the quantum state as there is no ponderomo¬ 
tive squeezing. At lower frequencies, rpond ^ r is huge 
and the injected vacuum is squeezed even more along the 
Os-axis and anti-squeezed along the Cc-axis. Then it is 
followed by a rotation by an angle ttpond = —7r/2 effec¬ 
tively making it a strongly amplitude squeezed vacuum, 
increasing quantum noise in phase quadrature by a factor 

g2(r-po„d+r) ^ I 

For Sagnac, the story is a bit more complex, as one has 
to account for non-zero rotation represented by Upond = 
Mpond -b 7r/2 and a much smaller ponderomotive squeez¬ 
ing factor, Tpond- Consider the case of no signal recycling 
with phase squeezed vacuum injection shown by dashed 
red lines in panels (b) and (f) of Fig. 6. First, it gets ro¬ 
tated counterclockwise by Upond — tt/S ~ 35° at LF and 
i^pond — 7r/4 ~ 45° at HF which is followed by a weak, 
J'pond — 3.25 dB ponderomotive squeezing (we look at 
low frequencies ~ 10 Hz). The latter stretches the ro¬ 
tated ellipse along the Uc-axis and squeezes it along the 


Os-axis by the same factor of 2. The resulting stretched 
ellipse inclines more towards the Og-axis (new angle 8°). 
Finally, it is rotated clockwise by Upond — —7r/3 ~ —55° 
at LF and Upond — —7r/4 ~ —45° at HF ending up 
with an ellipse at almost 45° to both axes, thereby hav¬ 
ing equal amount of noise in both, phase and ampli¬ 
tude quadrature. At higher frequencies quantum noise 
remains squeezed at phase quadrature. 

Two conclusions can be derived from the above qual¬ 
itative discussion: (i) because of very strong pondero¬ 
motive squeezing, Michelson interferometer can benefit 
from squeezing injection only if squeezing angle varies 
with frequency (frequency dependent squeezing), or the 
optimal readout quadrature is chosen at each frequency 
(variational readout), and (ii) in Sagnac interferometer, 
however, one can find the right values of squeezing fac¬ 
tor, squeezing angle and readout quadrature to partially 
compensate for ponderomotive squeezing. The resulting 
sensitivity curves for both configurations are shown in 
Fig. 6 and optimal parameters are given in Table H. 

For constant squeezing angle. A, and readout quadra¬ 
ture, (/, the optimal broadband configuration requires 
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FIG. 6: Quantum noise sensitivity of Michelson (left, panel (a)) and Sagnac (right, panel (b)) interferometers quantum in 
case of constant phase squeezed vacuum injection. Panels (c) and (d) show ponderomotive squeezing factor, rpond(fl), while 
panels (e) and (f) show quantum noise ellipse overall rotation angle, i.e. rot. angle = A + Mpond + Upond of Michelson and 
Sagnac interferometers, respectively. Solid black lines of quantum noise sensitivity of bare interferometers (w/o recycling and 
squeezing) are given for reference. Dashed lines demonstrate the effect of constant phase squeezing (r = 10 dB, A = 0) 
with phase quadrature readout = it/2), dash-dotted lines show quantum noise for optimised squeezing angle and readout 
quadrature. Parameters used for these plots are given in Table II. Light’s quantum state transformation by an interferometer 
with squeezed vacuum injection is illustrated by noise ellipses. Ellipses on panels (a) and (b) stand for injected squeezed vacuum 
noise ellipses (left one for 10 dB phase squeezed input, right one for optimised squeezed input). Ellipses on panels (e) and (f) 
demonstrate quantum state of the outgoing light transformed by ponderomotive squeezing at different frequencies (10, 100, 
and 1000 Hz). Green arrows point at the direction of readout quadrature, (. An upper row of ellipses on both lower panels 
refers to the 10 dB phase squeezed vacuum injection case. A lower row stands for the optimised fixed angle squeezed vacuum 
injection case. 


Parameter 

Notation 

Michelson 

Sagnac 

Mirror mass, kg 

m 

200 

200 

Arm length, km 

L 

10 

10 

Girc. power, MW 

Pc 

3.0 

3.0 

ITM transmittance, % 

Titm 

5 

9 

Squeezing factor, dB 

CdB 

5.6 

6.3 

Squeezing angle, deg. 

0sqz 

-6° 

15° 

Homodyne angle , deg. 

C 

83° 

97° 


TABLE II: Optimal parameters for configurations plotted in 
Fig. 6. 


squeezing factor to be limited. If A and C were optimal 
at each frequency, i. e. satisfied the following relation (cf. 


Eq. (16) of [48]) that delivers minimum to Eq. (39): 

cot^-I-tanA =/C(Q), (40) 

there would be no such limitation and the resulting quan¬ 
tum noise would decrease by the factor 


However, this can be achieved only by having either A, 
or C, or both frequency dependent to satisfy Eq. (40) at 
each frequency Q. 

b. Frequency dependent squeezing. Optimal fre¬ 
quency dependence can be followed quite closely using a 
phase dispersion of an ordinary optical cavity. This idea 
was raised by Kimble et al. [24] who suggested to in¬ 
troduce deliberately constructed optical cavities between 
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FIG. 7: Quantum noise of Michelson (left, panel (a)) and Sagnac (right, panel (b)) interferometers with frequency dependent 
squeezing injection: solid black lines show quantum noise without squeezing, thin dash-dotted lines show quantum noise 
with lossless single filter cavity with parameters given in Table III, and thick dashed lines show quantum noise with account 
for optical loss in the FC. Hatched regions demonstrate sensitivity deterioration due to FC loss. Panels (c) and (d) show 
ponderomotive squeezing factor vs. frequency, and panels (e) and (f) illustrate frequency dependent phase space rotation of the 
light noise ellipse in the interferometer at all stages, i.e. before the ponderomotive squeezing is applied (dashed lines), after it 
(dash-dotted lines), and the overall rotation angle (black dotted lines). Corresponding output noise ellipses are shown for both 
interferometers in panels (c) and (d) at 3 different frequencies (10, 100 and 1000 Hz). Readout quadrature direction, given 
by angle Q is illustrated by green arrows. Quantum noise at each frequency can be conveniently represented by a projection 
of the corresponding noise ellipse on the readout quadrature direction. These plots show that a single filter cavity is able to 
provide only a quasi-optimal phase rotation. The effect of this deviation from optimal dependence is most pronounced near the 
medium frequencies (around 100 Hz). Effect of optical loss in the arms is not negligible for Sagnac interferometer and shown 
by pink colour-filled area at low frequencies (Tbtm = 40 ppm). 


the dark port of the interferometer and a squeezed vac¬ 
uum generator, thus making frequency dependent A(Q), 
or before the readout stage to make frequency dependent 

m- 

The latter idea known as variational readout is more 
beneficial, in theory, as setting readout quadrature to 
Q — arccot[A(!(Q)] makes radiation pressure contribution 
to the final spectral density effectively zero (see, e.g., 
Eq. (29)). Closer look, accounting for all losses, has 
shown that this advantage vanishes and the modification 
of the alternative of frequency dependent squeezing in¬ 
jection has proven to pay off more (see Sec. 6.1 of review 
[31] and references therein). 

Omitting loss, filter cavity 1/0-relations can be written 
as: 

df{yt) = Tf{n)i = ( 42 ) 


where transfer matrix, T/(Q) = (Q), represents 

a simple phase shift, of no interest to us, plus phase 

plane rotation with frequency dependent rotation angle 

2 "^ ^1 (5 f 

0,(n) = arct.an-5— 

where 7/1 = cTf/{ALf) and Sf are the filter cavity half¬ 
bandwidth and detuning, respectively, and Tf stands for 
a filter input mirror power transmissivity while T/ is its 
length). 

Using (38), quantum state of the light after the fil¬ 
ter cavity can be written in the form of spectral density 
matrix as: 

S“ =K;,+,^,(Q)-S[2r].Kl+^/Q), (43) 

which can then be substituted into the general equation 
(17) to yield spectral density of the interferometer with 
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Parameter 

Notation Michelson Sagnac 

Mirror mass, kg 

m 

200 

200 

Arm length, km 

L 

10 

10 

Circ. power, MW 

Pc 

3.0 

3.0 

ITM transmittance, % 

Titm 

5 

9 

Squeezing factor, dB 

rdB 

10.0 

10.0 

Squeezing angle, deg. 

0sqz 

-11° 

-19° 

Homodyne angle , deg. 

C 

80° 

72° 

Filter cavity parameters 

FC detuning, Hz 

Sf 

57 

44 

rFC bandwidth, Hz 

7/1 

48 

117 

LFC mirr. transmittance, ppm/m 

Tf/Lf 

4 

10 

rFC loss bandwidth, Hz 

7/2 

12 

12 

LFC round trip loss, ppm/m 

^flLf 

1.00 

1.00 


TABLE III: Optimal parameters for configurations plotted in 
Fig. 7. 


frequency dependent squeezing injection. The result, in 
the absence of loss, will be a formula similar to (39) but 
with substitution A —>■ A + ^/. 


If filter cavity parameters, 7/1 and Sf are chosen prop¬ 
erly, as well as the initial squeezing angle. A, one gets 
a frequency dependent rotation of the injected noise el¬ 
lipse such as to almost completely compensate the pon- 
deromotive squeezing. Fig. 7 demonstrates how it hap¬ 
pens. Quantum noise spectral densities for Michelson 
and Sagnac interferometers with single input FC are 
shown (panels (a-b)). Transformation of input phase- 
squeezed vacuum state in the interferometers is illus¬ 
trated by noise ellipses for output light at different fre¬ 
quencies (10, 100 and 1000 Hz) and by frequency de¬ 
pendence of ponderomotive squeezing factors and state 
rotation angles. 


However, an overall quantum noise suppression shown 
in Fig. 7 is not homogeneous over frequencies. This is 
due to optical loss in the filter cavity, taken to be 1 
ppm/metre for these plots which is achievable by cur¬ 
rent technology [25]. Loss in the cavity causes an in¬ 
coherent mix in of vacuum photons pertaining thereto, 
with squeezed vacuum fields entering the cavity, thereby 
degrading its squeezing factor. This effect is most pro¬ 
nounced within the filter cavity bandwidth which ex¬ 
plains worse than expected improvement of quantum 
noise sensitivity at low frequencies (< 100 Hz) [28, 31, 
49]. The detailed account for optical loss in various ele¬ 
ments of the interferometer and its influence on the over¬ 
all sensitivity is briefly outlined in the section below and 
considered in detail in Appendix C 


G. Optical loss in gravitational wave 
interferometers 


As we argued in the Introduction, optical loss makes 
the major constraint on the achievable sub-SQL sensitiv¬ 
ity for any scheme of GW detector relying on quantum 
correlations for beating the SQL [50] , which includes use 
of squeezed light and optimal quadrature readout. Its in¬ 
fluence on detector sensitivity is quite diverse depending 
on where it originates in the optical layout. Nevertheless, 
the way quantum correlations are degraded by optical 
loss is the same. Incoherent vacuum fields that due to 
Fluctuation-Dissipation Theorem of Callen and Welton 
[30] accompany loss sources in any system admix with 
squeezed light travelling through the interferometer. 

c. Losses in the readout train. It is well understood 
that optical loss in the elements of readout optical train, 
including the non-unity quantum efficiency of photode¬ 
tectors, can be reduced to a single, frequency indepen¬ 
dent coefficient of equivalent quantum efficiency, rid = 
l — Cd < 1 , where < 1 can be thought of as a fractional 
photon loss at the photodetector [24, 28, 31]. Frequency 
dependence can be safely omitted here, for any resonant 
optical element in the readout train, including output 
mode cleaners (OMC), has bandwidth much larger than 
the detection band of the interferometer. 

The expression (16) for an output observable of the 
GW interferometer is modified in the presence of readout 
losses as follows: 


= v/l - Cd (be cos C + h sin c) 

+ (nd- c cos C + fid- s sin C) 

= -b+V^nJ-rid, (44) 

where fid = {fid- c, fid-, s}"*" stands quadrature vector of 
loss-associated vacuum fields with unity spectral density 
matrix. 

Spectral density formula (17) in lossy readout case will 
be simply: 


‘S'pDloss(^) — ^SQL 


nr 


[T-Sr 




(45) 


where ^d = \l<^dl{f - ^d)- 

d. Optical loss in a squeezing injection optics. Re¬ 
cent successes in generating low frequency strong squeez¬ 
ing and using it to suppress quantum noise of the GW 
interferometers GEO 600 and LIGO below the vacuum 
level [51, 52] have indicated optical loss in the injection 
train as the main barrier for getting highly squeezed vac¬ 
uum entering the GW detector dark port. A compli¬ 
cated, multi-element optical setup is required to make the 
squeezed light leaving a squeezer match perfectly with 
the spatial profile of the carrier mode of the interferome¬ 
ter. This results in reasonably high level of integral opti¬ 
cal loss in this installation, which can be characterised by 
an integral, frequency-independent injection power loss 













13 


coefficient, Cgqz- Following the same chain of argument 
as for the readout train loss, the 1/0-relation for the in¬ 
jection train can be written as: 


where i and d stand for input and output fields of the 
FC, respectively, and h represents vacuum fields due to 
loss. Transfer matrices for filter cavity are defined as: 


^dark port 




+ 


5 ns. 


(46) 


where idark port stands for the light field, entering the 
dark port of the detector (or the filter cavity in case of 
frequency dependent squeezing injection), and igqz and 
TZsqz are the field generated by a squeezer and a vacuum 
field due to injection loss, respectively. If the squeezer 
is capable of generating squeezed vacuum with e“’'(e’’) 
(anti-)squeezing factor, the effective (anti-)squeezing fac¬ 
tor at the dark port reads: 

g 2reft = + Csqz , 

^ (1 _ ( 47 ) 

e. Optical loss in the arms and in filter cavities. Op¬ 
tical loss in Fabry-Perot cavities, such as arm cavities 
and filter cavities, is known to have frequency depen¬ 
dence with the major impact at low sideband frequencies 
within the cavity optical bandwidth. A very illuminating 
discussion on this subject is given in [28] where optical 
loss in filter cavities is studied in detail. The main source 
of such loss in large suspended cavities is the scattering of 
light off the mirror surface imperfections of microscopic 
(micro-roughness) and relatively macroscopic (’’figure er¬ 
ror”) size [28, 53]. 

In general, optical loss in the cavity depends on the 
cavity length in an involved way (see Appendix C in [28]). 
However, for a fixed Hlter cavity length {Lf =40 metres 
in our case) and a fixed value of total photon loss per me¬ 
tre {ef = 1 ppm/m), a conventional description of Fabry- 
Perot cavity with one lossy mirror (usually, an ETM one) 
and another lossless one (an ITM, respectively), works 
perfectly Hne. A detailed derivation of lossy cavity I/O- 
relation is given in Appendix C 1. Here we only present 
its general form which reads: 


arm 

with Tarm = -|- a transfer matrix for input 

fields, darm, Narm = -|- Narm is & transfer matrix 

for loss-associated vacuum helds, narmj and tarm is an 
optomechanical response function of the cavity defined 
by Eq. (15). We wrote transfer matrices Tarm and Narm 
as sums of shot-noise component, and back- 

action component, T!]'j.^/Narm) to discern Fabry-Perot 
cavities with movable mirrors, as in the arms, from the 
ones with fixed mirrors, as in filter cavity. In the lat¬ 
ter case, the back-action components equal to zero, as 
well as the optomechanical response function tarm = 0. 
The exact expressions for these matrices are derived in 
Appendix Cl. In a special case of filter cavity, the I/O- 
relations can be simplified by omitting back-action and 
signal parts in the above formula: 

d/(r2) = T/i(H)-I-N/n, (49) 


= T 


arm'-*'arm 


-N, 


arm^arm 


+ ta 


dsQL 


(48) 


T; = (H), Nf = N|-(fl), (50) 

with expressions for and given by Eqs. (C4) 
and (C5), respectively. 

Using this simplihed formula, quantum noise spectral 
density for interferometer with lossy input filter cavities 
can be obtained by substituting into (45) the following 
expression for input field spectral density matrix: 

= T; • K;, • S[2r] • K; • T} + % • , (51) 

instead of (38). The last term here peaks near the res¬ 
onant frequency of the cavity which thereby decreases 
squeezing of the vacuum fields entering the cavity. But 
the off-resonant squeezed vacuum fields reflect off the FC 
without deterioration. This explains why optical loss in 
the cavities have major impact at low frequencies within 
the FC linewidth. The effect of lossy FC on interfer¬ 
ometer sensitivity is illustrated in panels (a) and (b) of 
Fig. 7 by hatched region between lossy and lossless FC 
sensitivity curves. 

/. Optical loss in Michelson and in Sagnac interfer¬ 
ometers Depending on the topology of an interferom¬ 
eter the influence of loss differs. In Michelson, the loss 
in each arm cavity sum up incoherently to produce ad¬ 
ditional, loss-induced quantum noise at the output port 
of the interferometer. The quantum noise I/O-relations 
for Michelson interferometer are the same to that of a 
single Fabry-Perot cavity,, a fact known as ’’scaling law” 
[31, 32]. Using I/O-relations at the beam splitter for the 
fields entering and leaving the interferometer at the dark 
port and for the fields at north (N) and east (E) arm 
cavities: 


p + i 


p — i 



(52) 


one can derive I/O-relations for lossy interferometer that 
read 


Omi = Tmi*- f NmiK + Imit-) (53) 

hSQL 

where 

Tmi = T arm 5 Nmi = N arm 7 Imi = arm 7 (54) 

and n = /y/2 is a loss-associated quantum 

noise. The only difference with a single Fabry-Perot is in 
a higher response to GW signal. 

Quantum noise spectral density matrix for a lossy 
Michelson interferometer can be written as 

= Tmi • s™ • Tj;,! + Nmi • N^j, (55) 
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and the resulting spectral density of the interferometer 
quantum noise in terms of GW metric variation reads: 


Qh ,loss 


m 


'^SQL 

IHWtMiP 


HT.{TMi-Sr 


(56) 


where we included losses in the readout train as well. The 
input state spectral density matrix, S“, is substituted 
either from Eq. (38) for hxed squeezing input, or from 
Eq. (51) for frequency dependent squeezing input with 
realistic filter cavity, or can be just an identity matrix 
for vacuum injection. 

In the case of small arm cavity round-trip loss, which 
can be reduced to ETM transmissivity Tetm ^ Titm, 
one can get the following closed expression for Michelson 
quantum noise spectral density: 




MI 


= S^. 


52 

'^SQL 




Af MI sin C 


+ CarmAv 


MI 


(57) 


where is the lossless quantum noise spectral den¬ 
sity defined in Eq. (39) with substitution K. JCmi, and 
^arm = VTetm /7?etm • Note that the last term in the ex¬ 
pression above is responsible for an additional radiation 
pressure created by the arm loss-associated helds n. This 
contribution is usually neglected, for it is much smaller 
than the main radiation pressure term in the Michelson, 
as — Tetm 10“® <C 1. In Sagnac, however, as we 
will see below, this approximation breaks down. 

In case of the Sagnac interferometer, two counter prop¬ 
agating light beams visit two arm cavities in sequence on 
their way to the output port of the interferometer. On 
each pass, a new loss-associated vacuum noise is gener¬ 
ated and contributes to both, shot noise and radiation 
pressure components of the total quantum noise. The 
I/O relations (25) shall be modified to: 

0 = Tsii -|- Ngjfii -f -I- tsiy , (58) 

hSQL 


where, for a Sagnac scheme based on ring arm cavities 

mi =-—, (59) 


ni = 


~LN ~RE 

n — n 


V2 


V2 


are the combined noise fields at the dark port originating 
from the first (I) and the second (II) pass. In case of 
polarisation-base Sagnac scheme, the above definitions of 
m and mi has to be multiplied by —1. Transfer matrices, 
Tsi, Ngj, Ng\ and signal response function, tgi are defined 
in Eqs. (C9), respectively. 

There is a signihcant difference between contributions 
to the sum quantum noise from the first pass, N|i, and 
from the second one, Ng\. This difference comes mainly 
from the different behaviour of back-action components 
of these two. For back action force created by the first- 
pass helds, ni, speedmeter radiation pressure suppres¬ 
sion, as in (6) in Sec. Ill A, works perfectly hne, as these 


helds share the same optical path as the dark port input 
held. For the second-pass helds, rin, this mechanism does 
not work as they only travel inside one arm cavity before 
leaving the interferometer through the dark port. This 
results in uncompensated Michelson-like radiation pres¬ 
sure force that shows up as a steeper (ex elevation 

of noise at low frequencies. 

Quantum noise spectral density matrix for a lossy 
Sagnac interferometer can be written as 

= Tsi • Sr • 4j + N^gj • (Nigi)t + • (N|\)t, (60) 

and the resulting spectral density of the interferometer 
quantum noise in terms of GW metric variation reads: 

= iH^tgiP ^ ■ {^SI • Sr • T^I+ 

+N^i-(Nigi)t+N^i .(Nig\)t}.Hc, (61) 

which can be written in a closed, simplihed form for the 
case of small loss (Aioss ^ Titm): 

(^h, loss _ Qh I ^SQL I i r2 iK 1 

-^SI - -^sqz ^-^ \ -jp -+ 4armNMI > , (62) 

2 [/Csism C J 

where is the lossless quantum noise spectral density 
dehned in Eq. (39) with substitution 1C —>■ /Csi- Note the 
same last term here, as in Eq. (57). Now this term can¬ 
not be neglected as /Cmi(^^ —)■ 0) ACsi(^ 0), which 
physically means that radiation pressure in Michelson is 
way too high compared to Sagnac interferometer, and 
even for fairly small arm loss, Tetm ~ 10“®, this term 
causes noticeable effect at low frequencies as shown by 
pink colour-filled area on the panel (b) of Fig. 7. 


H. Combined signal recycling and frequency 
dependent squeezing injection. 

Finally, we consider combined signal recycling and fre¬ 
quency dependent squeezing injection. Above we saw 
that Sagnac interferometer benefits a way more from ad¬ 
ditional filter cavity than Michelson, because of much 
weaker ponderomotive squeezing due to back-action. 
However, a non-perfect match of spectra of ponderomo¬ 
tive squeezing angles, Upond(f2) Mpond(f2), and that of 
the Hlter cavity phase transfer function, (j)f{Cl), results 
in worsened sensitivity at medium frequencies, Clj2TT ^ 
100 Hz. Panel (d) of Fig. 7 displays the result of this 
mismatch in a clear way as a non-optimal angle between 
the readout quadrature direction, /■, and the orientation 
of the noise ellipse at around this frequency. 

It turns out that for a Sagnac this mismatch can be 
reduced by means of detuned signal recycling. In this 
case, the non-monotonous dependence of ponderomotive 
angle on frequency allows for more effective use of the 
filter cavity (see panel ((d) in Fig. 8). Such a dependence 
comes from a structure of the detuned signal-recycled 
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FIG. 8: Quantum noise sensitivity of Michelson (left, panel (a)) and Sagnac (right, panel (b)) interferometers in case of 
combined frequency dependent squeezed vacuum injection and optimal signal recycling. Red vertical lines denote the location 
of two optical poles of the SI defined in Eq. (64). Solid black lines in panels (c) and (d) show ponderomotive squeezing in 
interferometers w/o signal recycling mirror, while red dashed lines demonstrate how it changes if an optimally detuned signal¬ 
recycling is employed. Panels (e) and (f) show frequency dependence of ponderomotive angles (red lines) and the optimal 
FC phase shift (black dashed curve) for a detuned signal-recycled interferometer. Light quantum state transformation by an 
interferometer with squeezed vacuum injection is illustrated by noise ellipses. Ellipses on panels (a) and (b) stand for injected 
squeezed vacuum noise ellipses for a 10 dB optimally squeezed vacuum. Ellipses on panels (c) and (d) demonstrate quantum 
state of the outgoing light transformed by ponderomotive squeezing at different frequencies (10, 100, and 1000 Hz). Green 
arrows point at the direction of readout quadrature, </■ Lower ellipse at 10 Hz represents a would-be output state, were there 
no optical loss in the filter cavity. All the relevant parameters for these plots are given in Table IV. 


Sagnac interferometer response function. As a system 
of two coupled Fabri-Perot cavities, it, quite expectedly, 
has two optical resonances manifesting themselves as two 
poles of both, the optical and the optomechanical transfer 
functions [54], that is {Tsi,sR; tsi,SR} oc I?“^(Q) where 


v{n) = (Q - Qi)(Q -b Q];)(Q - n2){^ + ni) , (63) 

f^l,2 = (5 i, 2 -b *71,2 , (64) 

2Vpsr cos Asr 


^1,2 — ^arm “b 4=7; 
7l,2 = 7arm 


1 -b PSR ± 27 ^ sin (/)SR 
1 - PSR 


1 -b PSR ± 2 ^psr sin (/)SR 


■710 


, (65) 

( 66 ) 


interferometers. However, unlike MI, the frequency de¬ 
pendent optical rigidity that arises in detuned Sagnac 
interferometer has no DC component and thus creates 
no mechanical pole in the mechanical response function. 
As shown in [54], it has the first non-vanishing compo¬ 
nent in the Taylor expansion that is oc thereby adding 
an inertial term to the test mass dynamics, i.e. 


-fbsi(H) = — 


m0armf2^(^l "b S 2 ) 

2V{n) 


^2 *R0arm(^l -b < 52 ) , 


and darm and 7arm are the arm cavities detuning and 
bandwidth and yioss = cAioss/(4L) is the additional 
bandwidth of the arm cavity due to a fractional round- 
trip photon loss Aarm therein. 

Similar to Michelson interferometer, the detuned signal 
recycling causes test mass dynamics variation in Sagnac 


where the coefficient in front of acts akin to additional 
mass, i.e. iFsi(^2) = —m.opt(D)Q^. This is the reason, 
why the detuned SI has no opto-mechanical pole (see 
Sec. 6.4.3 of [31] and references therein for more details). 

The prospects of using optical inertia in Sagnac re¬ 
mains an open question, for its sign changes between 
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Parameter 

Notation Michelson Sagnac 

Mirror mass, kg 

m 

200 

200 

Arm length, km 

L 

10 

10 

Circ. power, MW 

Pc 

3.0 

3.0 

ITM transmittance, % 

Titm 

5 

9 

Squeezing factor, dB 

rdB 

10.0 

10.0 

Squeezing angle, deg. 

^sqz 

-9° 

-36° 

Homodyne angle , deg. 

C 

81° 

55° 

Signal recycling cavity parameters 

SRM transmittance, % 

Tsrm 

12 

76 

SRC detuning, deg. 

<^SRC 

90° 

112° 


Filter cavity parameters 


FC detuning, Hz 

Sf 

13 

315 

rFC bandwidth, Hz 

7/1 

17 

307 

LFC mirr. transmittance, ppm/m 

Tf/Lf 

1.4 

26 

rFC loss bandwidth, Hz 

7/2 

12 

12 

LFC round trip loss, ppm/m 

^flLf 

1.00 

1.00 


TABLE IV: Optimal parameters for configurations plotted in 
Fig. 8. 



Frequency, Hz 


FIG. 9: Real part of the mechanical susceptibility of the 
Sagnac interferometer dARM mode with detuned signal re¬ 
cycling whose spectral density is plotted in the right panel 
of Fig. 8. Black dashed line represents the susceptibility of a 
free mass and, by coincidence, is the mechanical susceptibil¬ 
ity of the Michelson interferometer dARM. Red vertical lines 
denote the location of two optical poles of the SI. 


the optical resonances, jOil ^ O ^ l^ 2 |, which makes 
it hard to provide a stable broadband sensitivity im¬ 
provement. The broadband optimisation favours di¬ 
minishing optical inertia. In Fig. 9, we plot the (real 
part of) mechanical susceptibility of a detuned signal- 
recycled Sagnac interferometer, xsi(^^) = —(MdARM^^^ — 
along with the free mass susceptibility, 
Xf.m. = —pertaining to the RSE config¬ 
uration of the Michelson interferometer. Here, fJ-dARM = 
m/4 is the reduced mass of the dARM mode of the in¬ 
terferometer. The two curves diverge slightly only at 
frequencies below 10 Hz. This means the influence of 
the opto-mechanical interaction on the optical poles is 
negligible and the above expressions (64) can be used. 

Comparing panels (b) of Figs. 7 and 8, one can no¬ 
tice that the optimal location of the two optical poles 


Parameter 

Notation 

LIGO ET 

Low frequency detection limit, Hz 

/min 

5 1 

High frequency detection limit, kHz 

/max 

5 10 

Optical pump wavelength, nm 

Xp 

1064 

Mirrors mass, kg 

m 

40 200 

Interferometer arms length, km 

L 

4 10 

Circulating optical power, MW 

P 

arm 

s: 3 

Optical poser passes through BS, kW 

Pbs 

s: 80 

Squeezing, dB 

rdB 

s; 20 

Arm cavity round trip losses, ppm 

Tetm 

40 

Photodetector power losses, % 


1 

Squeezer power loss, % 

^sqz 

5 

FC round trip losses, ppm/m 

Af/Lf 

1, 0.1, 0.01 


TABLE V: The main parameters and their numerical values. 


encompasses the region of middle frequencies, where fre¬ 
quency dependent squeezing alone fails to compensate 
the ponderomotive squeezing rotation completely, yield¬ 
ing a hump-like behaviour of the quantum noise spectral 
density. Separation between optical poles provided by de¬ 
tuning allows to shape the frequency dependence of pon¬ 
deromotive rotation angles, Upond(^^) and Upond(i^), in 
such a way to match the filter cavity phase rotation char¬ 
acteristic almost perfectly. The result of this combined 
action is illustrated by the plots of Upond{^), 'Cpond(f^) 
and </>/(H), as well as noise ellipse diagrams in panel (d) 
of Fig. 8. One can see that the readout quadrature an¬ 
gle, now fits the orientation of the 100 Hz noise ellipse 
perfectly. 

Note also that for Michelson interferometer, the broad¬ 
band optimisation favours RSE configuration without ad¬ 
ditional rigidity and, therefore, without additional op¬ 
tomechanical pole. This can be explained by the un¬ 
welcome effect of the optical spring, turning the dARM 
mode into an oscillator, on the rotation of the noise el¬ 
lipse around the new mechanical resonance frequency it 
introduces. The ensuing frequency dependence Upond(f^) 
and Upond(^^) could hardly be matched by a single filter 
cavity. 


IV. NUMERICAL OPTIMISATION 
A. Optimisation protocol and figure of merit 

In this section we discuss the optimisation protocol 
we applied to obtain optimal broadband configurations 
described above. We start with the definition of the 
parameter space for each configuration. For a sim¬ 
ple signal-recycled interferometer w/o squeezing injec¬ 
tion the parameter space consists of four-parameter vec¬ 
tors PSR = {7Arm,RiTM,/0sR,</sR}, including arm cir¬ 
culating power, Farm, power reflectivity of the ITM, 
Fitm, and two parameters of the SR cavity, the SRM 
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amplitude reflectivity, psR the SRC phase shift, 
(pSR. Introduction of fixed input squeezing extends 
the parameter space by two more dimensions, namely 
squeezing factor, r, and angle. A, thus yielding Psqz = 
{-Parm,RiTM,PSR,</>SR,'r,A}. Finally, we include filter 
cavity to account for frequency-dependent squeezing in¬ 
jection, which adds two more parameters to the param¬ 
eter vector, textiti.e. the filter cavity detuning, 6f, and 
bandwidth, 7 /, giving the 8 -dimensional parameter vec¬ 
tor Pf = {Parm,7?ITM,PSR,</>SR,?', A,7 /i,5/}. 

Other parameters of the interferometer, such as arm 
length, mirror mass, loss coefficients etc. were fixed and 
their values used for optimisation are given in Table V. 
This table also outlines the constraints we put on the 
parameters of the interferometer. Another limitation we 
set is an upper bound of 80 kW on the power at the beam 
splitter, PbSj to mitigate detrimental effect of thermal 
lensing. As Pbs = 2Parm > this bound sets a 

limit on Pitm for given Pami and the other way round. 

The figure of merit we use for optimisation is the log¬ 
arithmic cost function for broadband detection proposed 
in [28], namely 


C(P) 



logio [' 5 ''‘( 27 r/;p) -f 5 ','[,f( 27 r/)] dlogio/- 


. . 

The set of parameters p that provides a global minimum 
for C'(p) is looked for numerically, using Nelder-Mead 
simplex method [55]. Here S'^(27r/;p) is quantum noise 
spectral density for the considered configuration of the in¬ 
terferometer and S^^f{2nf ) is a reference curve discussed 
above in Sec. II and defined in (5). By construction, this 
protocol seeks to minimise the area under the sum of 
these two curves in logarithmic scale. 


B. Results of optimisation. 

The main optimisation results are summarised in 
Fig. 1 and in Table I. We evaluated here the poten¬ 
tial of using Sagnac interferometer topology in conjunc¬ 
tion with other prospective quantum-noise-enhancement 
techniqnes, such as injection of frequency-dependent 
squeezed vacuum, detuned signal recycling and balanced 
homodyne detection. Results of our study clearly show 
that for all things being equal, Sagnac interferometer 
quantum noise has much higher potential for improve¬ 
ment compared to the Michelson topology. 

In this section, we focus on two particular questions: 
(i) how optimal sensitivity depends on the level of loss in 
the filter cavity and (ii) is there a significant benefit in 
changing mirror mass, thereby decreasing susceptibility 
of the interferometer to radiation pressure force. 

We vary loss of the filter cavity only for its impact on 
quantum noise is known to be dominating [25, 31, 49]. 
The results are presented in Figs. 10. Sagnac interferom¬ 
eter is proven by this result to be robust to FC loss in¬ 
fluence to much greater extent than its rival. The reason 


behind this result is quite simple. The adequate measure 
of optical loss in FC is the ratio of input mirror transmis¬ 
sivity to the round-trip absorption therein. Since much 
weaker back action in Sagnac interferometer allows for 
much broader band of FC than that that is required for 
Michelson interferometer, this ratio is much greater for 
the former one. Hence much weaker dependence on op¬ 
tical loss. Most clear illustration of this advantage is 
shown in Fig. 11 where the ratios of Michelson quantum 
noise amplitude spectral density and that of the Sagnac 
interferometer are plotted vs. GW frequency for different 
levels of optical loss. The larger is the loss the stronger is 
low-frequency improvement of the Sagnac interferometer 
over the Michelson one. Optimal parameters for consid¬ 
ered configurations are listed in Table VI. 

Variation of the mirror mass, to the contrary, does not 
favour any of the two considered configurations, as one 
can see in the Fig. 12. These plots clearly demonstrate 
that since the mechanical susceptibility of the interfer- 
ometr to an external force (including radiation pressure 
noise) changes identically for both, the Michelson and the 
Sagnac interferometers, the impact of mass variation on 
the overall quantum noise sensitivity must be the same 
for both of them. Optimal parameter values of consid¬ 
ered configurations are given in Table VH. 


V. CONCLUSION 

In our article, we conduct a detailed study of Sagnac 
interferometer as an alternative configuration for the 
3-rd generation GW detectors. Our analysis inclnded 
the most robust quantum non-demolition technologies 
planned to be used in the second generation detectors, 
namely, frequency-dependent squeezed vacuum injection 
and detuned signal recycling. We made a fair comparison 
of equivalent Michelson and Sagnac interferometers using 
basic parameters for two prospective projects of future 
detectors, namely, Advanced LIGO and Einstein Tele¬ 
scope. In both cases, we clearly demonstrated a much 
stronger potential for broadband quantum noise suppres¬ 
sion by Sagnac interferometer. We derived optimal pa¬ 
rameter sets for both projects. 

In case of aLIGO-like scheme, i.e. for arm length of 
4 km and maximum mirror mass of 40 kg, the quantum 
noise of Sagnac interferometer is consistently better by 
factor of 3 to 10 times than the quantum noise of the cur¬ 
rently planned baseline aLIGO broadband configuration 
(see solid red line in left panels (a) and (c) of Fig. 1). 
Michelson interferometer with the same base parameters 
shows much poorer performance at lower frequencies, as 
shown by blue dashed lines in left panels (a) and (c) of 
Fig. 1. 

For a more ambitious European project Einstein 
Telescope that features xylophone configuration of two 
Michelson interferometers which individual sensitivities 
are tailored so, as to provide best low- and high-frequency 
performance, respectively, we demonstrate that a single 
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FIG. 10: Quantum noise spectral densities for LIGO-like (upper panels) and ET-like (lower panels) interferometers with different 
levels of optical loss in the filter cavity. Parameters for all shown configurations are listed in Table VI. 
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FIG. 11: Improvement factor of the Sagnac interferometer with respect to the Michelson interferometer for different hlter cavity 
loss levels. Left panel: ET-like detectors (cf. lower panels of Fig. 10); Right panel: LIGO-like detector (cf. upper panels of 
Fig. 10). 


Sagnac interferometer can do almost equally well, or even 
better in terms of quantum noise for much lower price. 
This is clearly seen from the solid red lines in the right 
panels (b) and (d) of Fig. 1. 

Two interesting conclusions can be derived from the 
comparison of optimal parameters for Michelson and 


Sagnac interferometers given in Table I. Firstly, due to in¬ 
herently low radiation pressure noise at low frequencies, 
requirements to filter cavities for Sagnac configuration 
are much relaxed compared to that for Michelson one. In 
the former case the optimal filter cavity has more than 
an order of magnitude broader bandwidth (thus lower 
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FIG. 12: Quantum noise spectral densities for LIGO-like (upper panels) and ET-like (lower panels) interferometers with different 
mirror masses. Parameters for all shown configurations are listed in Table VII. 


finesse) than in the latter one, which results in much 
weaker dependence of sensitivity on optical loss in these 
devices (see Fig. 10). Secondly, strong radiation pressure 
at low frequencies precludes Michelson interferometer to 
benefit from stronger squeezing (optimal squeezing fac¬ 
tor is lower than the upper limit allowed by the optimi¬ 
sation). to the contrary, Sagnac interferometer requires 
as strong squeezing as possible, thereby making full use 
of this experimentally most developed QND technology. 


Finally, we have shown that slightly detuned from 
resonance, low-finesse signal recycling cavity allows to 
improve medium-frequency quantum-noise sensitivity of 
Sagnac interferometer by tailoring the ponderomotive ro¬ 
tation angles frequency dependence so as to match the 
filter cavity’s phase response function, thereby providing 
almost optimal output squeezed state rotation angle. 
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Appendix A: SQL for multi-mirror mechanical 
modes of an interferometer 

Consider an effective mechanical mode of the inter¬ 
ferometer that can be defined in most common form as: 

= (Ai) 

j=l 1 
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Parameter 

Notation 

LIGO 

Michelson 

Sagnac 


Michelson 

ET 

Sagnac 


Mirror mass, kg 

m 

40 

40 

40 

40 

40 

40 

50 

100 

300 

50 

100 

300 

Arm length, km 

L 

4 

4 

4 

4 

4 

4 

10 

10 

10 

10 

10 

10 

Circ. power, MW 

Pc 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

3.0 

3.0 

3.0 

1.8 

2.2 

2.8 

ITM transmittance, % 

Titm 

15 

15 

15 

8 

7 

7 

5 

5 

5 

16 

14 

11 

Squeezing factor, dB 

r-dB 

13.2 

19.4 

20.0 

17.9 

18.6 

18.7 

13.3 

13.2 

13.5 

18.5 

18.2 

18.2 

Squeezing angle, deg. 

0sqz 

-6° 

-9° 

-10° 

-17° 

-12° 

-11° 

-7° 

-5° 

-4° 

-13° 

-13° 

-11° 

Homodyne angle , deg. 

C 

O 

00 

81° 

00 

O 

o 

O 

CO 

00 

o 

79° 

O 

CO 

00 

O 

uo 

00 

86° 

O 

00 

77° 

79° 


Signal recycling cavity parameters 


SRM transmittance, % 

Tsrm 

82 82 70 

89 95 96 

9 11 13 

93 90 

87 

SRC detuning, deg. 

<^SRC 

90° 90° 90° 

101° 104° 105° 

90° 90° 90° 

100° 99° 

98° 


Filter cavity parameters 


FC detuning, Hz 


34 

33 

28 

269 

172 

154 

24 

18 

10 

155 

157 

147 

FC bandwidth, Hz 

7/1 

34 

28 

24 

806 

587 

544 

25 

21 

16 

644 

600 

539 

-FC mirr. transmittance, ppm/m 

Tf/Lf 

2.8 

2.3 

2.0 

68 

49 

46 

2.1 

1.7 

1.3 

54 

50 

45 

FC loss bandwidth, Hz 

7/2 

12 

1.2 

0.1 

12 

1.2 

0.1 

12 

12 

12 

12 

12 

12 

-FC round trip loss, ppm/m 

Af/L; 

1.00 

0.10 

0.01 

1.00 

0.10 

0.01 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


TABLE VI: Optimal parameters for LIGO-like and ET-like detectors for different values of optical loss in filter cavity. 


Parameter 

Notation 

LIGO 

Michelson Sagnac 


Michelson 


ET 

Sagnac 


Mirror mass, kg 

m 

100 

200 

100 

200 

50 

100 

300 

50 

100 

300 

Arm length, km 

L 

4 

4 

4 

4 

10 

10 

10 

10 

10 

10 

Circ. power, MW 

Pc 

1.0 

1.0 

1.0 

1.0 

3.0 

3.0 

3.0 

1.8 

2.2 

2.8 

ITM transmittance, % 

Titm 

15 

15 

6 

5 

5 

5 

5 

16 

14 

11 

Squeezing factor, dB 

rdB 

13.0 

13.1 

18.1 

18.3 

13.3 

13.2 

13.5 

18.5 

18.2 

18.2 

Squeezing angle, deg. 

0sqz 

-4° 

-3° 

-15° 

-13° 

-7° 

-5° 

-4° 

-13° 

-13° 

-11° 

Homodyne angle , deg. 

C 

O 

00 

00 

o 

75° 

77° 

O 

CO 

00 

85° 

86° 

O 

00 

77° 

79° 


Signal recycling cavity parameters 


SRM transmittance, % 

Tsrm 

88 

91 

87 

85 

9 11 13 

93 

90 

87 

SRC detuning, deg. 

4>sb.c 

90° 

90° 

99° 

98° 

90° 90° 90° 

100° 

99° 

98° 


Filter cavity parameters 


FC detuning, Hz 

<5/ 

23 

17 

217 

179 

24 

18 

10 

155 

157 

147 

FC bandwidth, Hz 

7/1 

26 

22 

710 

648 

25 

21 

16 

644 

600 

539 

-FC mirr. transmittance, ppm/m 

Tf/Lf 

2.2 

1.9 

60 

54 

2.1 

1.7 

1.3 

54 

50 

45 

FC loss bandwidth, Hz 

7/2 

12 

12 

12 

12 

12 

12 

12 

12 

12 

12 

-FC round trip loss, ppm/m 

Af/Lf 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


TABLE VII: Optimal parameters for LIGO-like and ET-like detectors with different mirrors’ masses. 


with Xj standing for individual displacements of each of 
the N interferometer’s moving mirrors with individual 
masses Mj. Coefficients aj are arbitrary constants dif¬ 
ferent authors prefer to choose individually (in most LSC 
papers it is assumed 1). Note that all Xj in this formula 
are taken with sign unlike the definition of a^dARM 
in the footnote on Page 3. This is simply because we 
choose a;-axis for each mirror so, as to absorb the sign 
information in the definition of Xj. 


Assuming all test masses free, one writes their equa¬ 
tions of motion in frequency domain as x“g(Il)x(n) = 
Fj with Xj-,o{^) = is the mechanical suscep¬ 

tibility of a free mass and Fj is the external force acting 
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on it. The equation of motion for then reads: 

where Xf^-,o{^) = > 


3 = 1 
N 


(A2) 


Xko{^)^^ 0-3 


Here index 0 in Xp.-, o indicates that this is the mechani¬ 
cal susceptibility of the mode x^, without any rigidities, 
either mechanical, or optical. By definition “ 

—thus yielding the following expression for the re¬ 
duced mass /i 


M = 




-1 


(A3) 


In GW interferometers, when only quantam noise is 
considered, external force Fj comprises of two compo¬ 
nents, the GW-induced one 

G,(n) = x-],(n)-^^^, (A4) 

with h the GW-strain amplitude, and the radiation pres¬ 
sure one (or back-action one) that can be written as: 

= Ff^- «(H) + K,{n) %(H), 

where F^-^- ^ is a fluctuational part of back action force 
and Kj is an optical rigidity. The requirement that all 
Kj shall be combined into a single optical rigidity of 
the effective mode sets a constraint on coefficients aj: 

Qj _ Xj; o(^) ^3 (^) w- ■ 

a,~ x^A^)K^[n)^ 

This is why, for instance, for a Fabry-Perot arm cavities 
with high-reflectivity mirrors all ai are equal, i.e. o.i = a. 

It can be shown that expressions (Al)-(A3) conserve 
the form of the system Lagrangian when it is rewritten 
for the effective mode. The conservation of traditional 
form of the Robertson-Schrodinger uncertainty relation 
[see eq. (148) in [31]] is one more indication of this fact: 

5,,(H)^Ff(H) - |5,F(f^)|% y , 

Here Sxx^ Sff and SxF are single-sided spectral den¬ 
sities of displacement measurement noise, back action 
noise and their cross-correlation spectral density for the 
effective mode x^, correspondingly. Therefore one can 
say that this single degree of freedom describes collective 
motion of all the test masses correctly. 

The single-sided force SQL spectral density for 
reads: 


where susceptibility function Xm(^) = 

[x-Vf^) + iX^(f^)]”' now includes optical rigidity 
Ar^(H). To get the SQL spectral density normalised to 
GW-strain, S'qql, one must calculate the factor .S, which 
relates tidal force G^(H) to the GW-strain amplitude 

h{n): 


G^(H) =il(H)h(H), 


’^sql(^) 


*^sql(^) 


In this paper, we assume aj = a, each mirror mass 
Mj = M, LjIemt = L, where L stands for arms length, 
and GjjiMT = 0. Hence, 

i? 2 xa = -^MLH 2 , 


which follows from eq. (A2), (A4). 

In case of free system, with XtJ.{^) = XtJ.-,o{^), 

and stands the free mass M force SQL by /|ql = 
'^fm sql(^) following: 


-^rxM SQl(^) 



/sQL ~ . 


Therefore, the GW-strain single-sided SQL spectral den¬ 
sity reads: 

4 

Qh /O'! _ ■'SQL _ ,2 

‘^4xMSQlI“1 — ~ ""SgL 


where hsQL is defined in Eq. (2). Note that GW-strain 
SQL does not depend on exact mechanical mode selec¬ 
tion, i.e. on a. 


Appendix B: Ponderomotive squeezing in GW 
interferometers 


Ponderomotive squeezing that takes place in a tuned 
lossless Michelson interferometer can be written as a se¬ 
quence of 3 unitary transformations - rotation, squeezing 
and second rotation [24]: 

[out) =e ^R(4ipond)*^(^pond)-^(^pond) ■ (® 1 ) 

where /3 is a scheme-specihc complex frequency- 
dependent phase shift which does not change the noise 
spectral density, the rotation operator R{a) and the 
squeezing operator S{r) are defined in Section 3.2 of 
[31]. Action of these operators on the vector of light 
quadratures, a = {fii, 02 }"'", results in a new vector, 
b = {bi, 62 }"'", that reads: 

b = Ta = K[Mpond] §[rpond] R^pond] a , (B2) 

with K the rotation matrix dehned in Eq. (32) and S the 
squeezing matrix equal to: 


e’' 0 

0 e-’' 


si<^^{n) = 2h\xFm-^ 


S[r] = 


(B3) 
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For a general optomechanical system without loss, the 
transfer matrix (TM) has a specific structure, namely, 
the optical TM is T™®®"® = and the radiation 

pressure one in proportional to oc t (crit)^, where a\ 
is the Pauli's matrix, also known as Ox ■ This structure of 
TM guarantees that the covariance matrix of the output 
field, V{, = TVqTI, will remain real like the input Vq, 
which are corresponded to gaussian \in) and \out) states. 
Factoring out common complex phase one ends up 
with a real matrix T^® = j']pmeas_|_']pb.a.j ^ singular 
value decomposition of which can be written as: 

T^® = K[Upond] S[T-pond] R^pond] , 

that proves that Eq. (B2) is indeed correct. 

In order to get the expressions for rpondi Wpond and 
■ypond, one can expand T^® in Pauli matrices: 

T^® = zqI + Ziai + Z2(72 + 


where 2 o,i, 2,3 are complex coefficients. 

Symmetries of the TM immediately allow to see that 
Z 3 = 0 and the zq = T|^^® = Tf°. Since all elements of T^® 
are real, the following relations hold for the remaining 
coefficients: 


opRe _|_ opRe opRe _ opRe 

cs sc ^ • cs sc 

Zi = --- = Zi, Z2 = I -T- =1' Z2, 


which means Zi, Z 2 are real. 

Then singular values can be calculated: 


Sl,2 


kl I ± 




Assuming = max{si,S 2 } and = 

min{si,S 2 } (i.e. rpond > 0 ) one can get the following 
expression: 


sinh Tpond 




if detT^® = 1 , 


z^ + zi, ifdetT^® = -l. 


1. Fabry-Perot interferometer with moving mirror 

I/O relations for a single Fabry-Perot arm cavity with¬ 
out any additional assumptions about its bandwidth can 
be written as: 


6a,rm = T, 


arm*-*'arm 


■N; 


arm^arm 


V^/iSQL ’ 


(Cl) 


where transfer matrices and signal response function 
read: 


Tarm — V^ITM^arm lPr„„'irETM]Pr„„NiTM - a/R ITm2 5 

^^arni = -v/TitmM arm lPr„„NETM , 
farm = \/2 TlTIVlMarmlPT^^tETM , 


Marm = P “ Pra™TETMPT„„TiTM] ^ , 

with 


(C2) 


Tetm 

Netm 

Metm 

Ietm 


\/Retm (I + Metm) , 
\/Tetm (I + Metm) , 


0 


0 


—RetmA^tm 0 
•\/2RetmA^tm 


(C3) 


Here = e*^'^“™M[cdprarm] describes the free propa¬ 
gation of light between the mirrors of the arm cavity. 


SiOpPc 
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arm'^arm 


is an optomechanical coupling factor for a single perfectly 
reflective free mirror and Pc stands for the full power of 
incident light beam. Note that here /if.p. sql = V^hsQL 
stands for the SQL of a Fabry-Perot cavity with 2 mov¬ 
able mirrors of mass M each, thus the factor \J2 in front. 


The expression for angles rtpond and Tpond are: 


a. Tuned cavity 


Mpond = arctan — - sgn [zi] y , 
Z Zq 4 

1 Z 2 r 1 ’’’ 

iipond = -- arctan-h sgn [zi\ - . 

Z Zq 4 


Appendix C: Transfer matrices for interferometers 
with losses 


Here we derive the I/O-relations for considered inter¬ 
ferometers beyond narrow-band approximation. 


In the case of resonance tuned cavity, the above ex¬ 
pressions simplify significantly and can be written as: 


T — 

arm — ^ ^ 




^arm — 


1 0 
—a/RetmA^^ arm 1 

1 0 
-M 1 

4 Re,tmK- am 
1 -I- Ritm 
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where in small optical losses approximation (i?ETM 1): 
(1 + -Ritm) TitmAI^tm 


^arm — 


1 — 2-y/ i?ITM COS 2S2r + i?ITM 


= arctan 


1 


tanflr 


N = 


N = 


T = t = l, 


1 — V i?lTM 

/1 — i?ETM ^ arm 

1 + i?ITM ^TM ’ 

/ Al^armAl^TM^ETM 1 + 

1 _ f?2_ p-i/3arm+inTarm 

r ^ITM c 


Following the reasoning after Eq. 48 one can derive 
the following expressions for shot-noise and back-action 
components of the transfer matrices: 


T” 


_ g2i/3arm][ 'jj'b.a. _ g2i/3a 


0 

—a/Set^/C 

a 


|i^meas _ g2/3armjyj^ p^b.a. _ gj/3armjy- 


0 0 
-Af 0 


0 
0 

(C4) 

(C5) 


b. Detuned narrow-band filter cavity. 


In case of filter cavities, mirrors can be assumed fixed 
and no radiation pressure effects are to be considered due 
to an absence of any significant classical light component 
therein. One can also make a so called narrow-band ap¬ 
proximation, assuming VlLf/c <C 1, where E/ is the filter 
cavity length and T/ <C 1 is input mirror power trans¬ 
missivity. Then one can write transfer matrices as: 


Tfc = 


Nfc = 


—h 


= 7/1 - 7/2 - (5/ + -h 2i07/2, 
h = —27/1(5/, 


2v'7/i7/2 


V 


7 / — ifli —6f 
Sf jf — in 


= 7/1 + 7/2 is a full 


where V = ( 7 / — Sj, 7 / 

cavity half-bandwidth and (5/ if its detuning. Here 7/1 = 
cTf/iALf) is a half-bandwidth part depending on input 
mirror transmissivity and 7/2 = c-Af/^ALf) is the loss- 
associated part of bandwidth with H/ <C 1 being the 
total round-trip fractional photon loss. 


2. Michelson/Fabry-Perot interferometer 

a. Fabry-Perot-Michelson (FPM) interferometer w/o 
signal recycling. 

I/O-relations of a Michelson/Fabry-Perot interferome¬ 
ter can be obtained by completing the above ones for the 
single arm with junction relations at the beam splitter: 

» ~ -N ~E 

-N i> + i -E P-i - -b 


E ■^N E •'N>E •'N,E 

where a ’ = a^rni) b = b^^^ stand for the input and 

output fields of the N and E arms, respectively. Hence, 
the Michelson interferometer I/O-relations read: 

o = TMli + NmiK + tMlT- ; (C6) 

hSQL 


where 


— Tfarm : I^MI — I^arni 7 ^MI — farm ■ 

Here n = ^narm — ^arm) /represents effective vacuum 
fields associated with optical loss in the arm cavities. 

In case of small losses the interferometer is described 
by opto-mechanical factor JCmi = l^arm and phase /3mi = 

/5arm ■ 


b. Signal-recycled Fabry-Perot-Michelson (FPM) 
interferometer. 

Along the same lines as is done in Sec. HID, the I/O- 
relations of signal recycling interferometer can be ob¬ 
tained from the following equations written for light fields 
on a signal recycling mirror (SRM): 


OSR = 


1 “ /^Sr'^SRO + PSRJPaisR 


* = PsR PSrPsrO + V 1 - PsR'P’a*SR 1 , 


(07) 


where an additional phase shift asR is introduced to 
satisfy the Scaling Law of [32] which maps the signal- 
recycled FPM interferometer and a single detuned Fabry- 
Perot cavity: 


PsR = e*- 


</sR = 


Wp^SR 


~ K[asR], QfSR = arctan 


PSR - 1 


tan 4>sn . 


, PSR + 1 

The solution of (06) and (07) gives the following: 

h 


osR = Tmi sr/sr + NmI SR« + ImI SR 


^SQL ’ 


where 


Tmi SR = Pet [(I — Psr)PsrMmisr'PmiPsr 
- PsrI] Pa , 

NmISR = \/l — PsRPaPsRMMISRN arm 5 

ImISR = \/l — PsRPaPsRMMISRtMI , 
MmISR = P — PSRTarmPsR] 
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3. Sagnac interferometer 

Full I/O-relation for the Sagnac interferometer are 
given in (58) and we repeat them here: 

o = Tsii + Ngifii + Ng’inii + tsiT-^ ■ 

hSQL 


b. Ring-arm-cavities Sagnac interferometer with signal 
recycling. 

To get now the transfer matrices and response for a 
signal recycled Sagnac interferometer, one have to solve 
equations (C7) and (58) jointly. This leads to the follow¬ 
ing expressions: 


The I/O-relations for each arm cavity of the Sagnac in¬ 
terferometer in general case read: 

i> = Tarm® + Narm® +(~ 1 ) tarm T- 1 " 

hSQL 

+ (C8) 

which accounts for an additional back action introduced 
by a counter-propagating beam into the response matrix 
of the beam under study. Here (—1)^ = 1, (—1)^ = 
— 1 and = narai- The effective vacuum fields that 
correspond to the optical losses in the arms in case of 
ring-cavity Sagnac interferometer are given by Eq. (59). 
For the polarisation Sagnac interferometer the definition 
of rii and fin differs only by the sign (minus). 


TsISR = Pa [(1 — Psr)PsrMsiSr'PsiPsR 

- PSRljPa , 


(CIO) 

PsR'PaPsRMsiSRNg’l” ’ 

(Cll) 

tsiSR = -(/l - PsR'^a'^SRMsiSRtsi ) 

(C12) 

where MgisR = P - PSRPaTsiPiR,] . 

(C13) 


Here again, Pq, ~ IR[agR], where ogR = 
arctan tanc^igR^ is the auxiliary phase shift 

added to satisfy the Scaling Law theorem [32]. 


4. Polarisation Sagnac interferometer with 
polarisation beam splitter (PBS) leakage 


a. Ring-arm-cavities Sagnac interferometer w/o signal 
recycling. 


Starting with the case w/o signal recycling, one can 
calculate transfer matrices using standard techniques and 
the result reads: 

TSI = — TarmMsiTarm , 

N^gl = Na™ - Ta,a,MgiN^,(;; , 

Ng\ = - TarmMgiNarm , (C9) 

tgi = —\/2 (I — TarmMgi) tarm , 

Mgi= [l + T^-]-'=I-Tb-. 

For small optical loss, i.e. i?ETM ^ 1, the SI transfer 
matrix and response function read: 


Tgi = 


1 0 
—-\/i?ETM^SI 0 


tsi = e 


_ piPsi 


4i?ETMAl^SI 
1 -I- i?iTM 


In this Appendix, we analyse in detail the input out¬ 
put relations of a polarisation Sagnac interferometer with 
lossy PBS, following the analysis of Wang et al. [36], 
but expanding their analysis to arbitrary signal-recycled 
Sagnac interferometers not limited by the approximation 
of narrow band of arm cavities. 

An ideal PBS is an optical element transmits 100% of 
p-polarised light and reflects by 90° the orthogonally s- 
polarised one. However, in reality, some fraction of the 
s-component gets transmitted and some fraction of the 
p-component gets reflected. To characterise this leak¬ 
age one usually introduces leakage coefficients, for 
reflected fraction of the p-component, and for trans¬ 
mitted fraction of the s-component. Physically, Cp and 
show the percentage of unwanted polarised light power 
mixed into the desired one. Keeping to the system of 
notations for the clockwise and counterclockwise propa¬ 
gating beams adopted in the main text, one denotes helds 
leaving the beam splitter toward the main interferometer 

^R,L 

and PBS as a , and the ones leaving it and mixing at 
the beam splitter to finally leave the interferometer for 

rJi,L 

the readout port as h . Hence we can get the following 
junction equations at the BS: 


where /3si and /Csi are given by Eqs. (27) and (28). Note 
the difference with the Michelson case, i.e. the double 
amount of power circulating in the arms of the Sagnac 
interferometer compred to that of an equivalent Michel¬ 
son one. This results in an additional factor of 2 in the 
definition of the circulating power in the interferometer. 
Pc = 4Parm; where Harm IS Understood as the single op¬ 
tical beam power, circulating in the arm cavity. 


^R 

_ Pp,s T —f*- _ Pp,s ~ _ ^p,s ~ ^p,s 

" V2 ’ " V2 ’ “ V2 

where for / J = {RE, LN} beams one has 


IJ /- 1 ^ j --^iJ 

d = VI - Csd +V^s> bp = y/l-epb -\-y/^p, 
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-LE 

b 


FIG. 13: Schematics of the held transformations on a lossy 
beam splitter (BS) and polarisation beam splitter (PBS). The 
straight arrows correspond to total helds of orthogonal polar¬ 
isations, while the arrows with semi-circular tails represent 
the two components comprising the outgoing held, i.e. the 
nominal polarisation component and the leaked component. 


and for IJ = {LE, RN} 


a. Polarisation Sagnac interferometer w/o signal recycling. 

Transfer matrices and response functions for both po¬ 
larisations can be shown to have the following form: 

= (1 - e,) Ml - yr^MsiT^^] - , 

= e,Mi - ycpMsiTL™] - , 

Til = Tgf = - epv^MiMsiTarm , 

Ni? = ^r^Mi [MgjNarm " , 

Nil = y/T^Mi - vT^MsiTarmNarm] , 

= yilMi [M^iNarm " y^armN^,^] , 

= v^Mi - v^MsiTarmNarm] , 

t|j = -^2 (1 - ep)Mi [I - \/l - EsTarm] tarm , 
tgl = -y^e^Mi [I - yi^Tarm] tarm ■ 

where 

Ml = [I - a/ 1 - Csy^MsiT^rm] > 

Msi = I - (v^ + VT^s) , 

The following relations that follow from the specific struc¬ 
ture of transfer matrices can be used to simplify the final 
expressions: 

rrpb.a. rrpb.a. ^b.a.i^ab.a. pipb.a. ii^b.a. ii^b.a. pipb.a. r\ 

armarm armarm armarm -^^arm-^arm ^ ? 

M|jM^i = M^iMPj = Msi, 

[Tarm, = [Mgi,Tarm] = [Mgi,Narm] = 0, 

MiiTb-=T^rm: , M^jtarm = tarm 


'IJ 


\/i 


—rrl rIJ 

Epflp T , 


bs 




e.,a„ 




IJ 


where Mgj stands either for Mgj, or Mgj, or Msi. 


The above expressions together with I/O-relations for the 
arms (08) give us the full system of equations whose 
solution can be written in the following form with P = 
{s,p}-. 


b. Signal-recycled polarisation Sagnac interferometer. 


Op — Tgf ip -|- Tgf ip Ngj^ hp -\- Ngj^np -I- tsiy- . 

h-SQL 

The correspondence of the above expression to the loss¬ 
less Sagnac interferometer I/O-relations of Eq. (58) 
is straightforward. The “signal” polarisation is p- 
polarisation, therefore ip and bp serve as an input and an 
output fields of the interferometer, and hp^ ^ corresponds 
to npii. 


Now we can also write the exact expressions for the 
transfer matrices and response functions of the signal- 
recycled Sagnac interferometer with losses and imperfect 
PBS: 


OSR. P — Tgi^RisR. P + '^Sl^SR^SR, P + 


PP 


Nsi^r^p + Ngi^j^rip -I- tsisR 


PP 


^SQL ’ 

































26 


where 


= P4(l - Psr) PsrC^^ - psrI ], 
"^SI^R = (l ~ ^'Sr) PoPsrC’^^Pq , 


N 


PR _ 
SI SR — 


\/l ~ PSR^^a'PsRMp 



+ PSRT^^'PsRMpNf« , 


t 


p _ 
SI SR — 


~ Psr'P“'PsrMp 




+ PSR'P'^^PsRMptgj , 


f^PR ^ 

Mp = 

f = 


PrmSRTr PsrT 

n -1 


TPR 

I - psrT^^PsrT^^Psr 
MpT^^PsR, 


PR 


Mp = [I - psrT^^PIr] 


and {P, R} = {p, s} in all combinations. 
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